NASA Technical Memorandum 4721 


Least Squares Best Fit Method for the Three 
Parameter Weibull Distribution: Analysis 
of Tensile and Bend Specimens With 
Volume or Surface Flaw Failure 


Bernard Gross 


JULY 1996 



National Aeronautics and 
Space Administration 


NASA Technical Memorandum 4721 


Least Squares Best Fit Method for the Three 
Parameter Weibull Distribution: Analysis 
of Tensile and Bend Specimens With 
Volume or Surface Flaw Failure 


Bernard Gross 
Lewis Research Center 
Cleveland, Ohio 



National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Program 

1996 


Least Squares Best Fit Method for the Three Parameter Weibull Distribution: 
Analysis of Tensile and Bend Specimens With Volume or Surface Flaw Failure 

Bernard Gross 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Summary 

Material characterization parameters obtained from natu- 
rally flawed specimens are necessary for reliability evalu- 
ation of nondeterministic advanced ceramic structural 
components. The least squares best fit method is applied to 
the three parameter uniaxial Weibull model to obtain the 
material parameters from experimental tests on volume or 
surface flawed specimens subjected to pure tension, pure 
bending, four point or three point loading. Several illustrative 
example problems are provided. 


Introduction 

The objective of this report is to apply the least squares 
best fit (LSBF) method to evaluate the parameters used in the 
uniaxial Weibull three parameter model. These parameters, 
scale factor a D , Weibull modulus m, and threshold (location) 
parameter are material dependent. Weibull two or three 
parameter models are used to specify a probabilistic distribu- 
tion for monolithic ceramic materials. The success in the use 
of the two parameter model rather than the three parameter 
model depends on the importance of ignoring the threshold 
(location) parameter. Disregarding the threshold parameter is 
conservative and simplifies matters. This simplification can 
be justified only by comparing the predicted behavior of a 
component with its observed performance. 

Equations are developed to obtain the three material 
parameters from inert volume or surface flawed data. Inert 
data imply fast fracture (no subcritical crack growth). The 
inert data are obtained from experimental tests on specimens 
subjected to either pure tension, pure bending, and four or 
three point loading (fig. 1). Ideally the data are obtained 
under conditions representative of the service environment. 

Several applications are presented in the section entitled 
EXPERIMENTAL APPLICATIONS. Experimental data are 
analyzed for volume flaw failure of silicon nitride (SNW- 
1000) specimens tested in four point bending (ref. 1). In 
addition, analysis is made of surface flaw failure data of 
silicon carbide specimens, annealed in both the longitudinal 
and transverse direction and tested in three point bending 
(Private communication from Sung Choi and Jonathan 
Salem, NASA Lewis Research Center). The four point bend 
volume flaw data are also used for a four point bend surface 
flaw analysis to illustrate the application of the developed 


equations. It is realized that these data are not representative 
of the physical problem. 

Symbols 

A tensile surface area 

b specimen thickness 

L beam length 

m Weibull modulus (shape parameter) 

n number of inert data points 

Pf probability of failure 

V volume in tension 

W specimen depth 

x,y,z Cartesian coordinates 

l° wer limit of integral for specimen 

(x,y,z) stress distribution in the specimen j at fracture 

0f Jmax maximum principal stress in specimen j at 

fracture 

°fJmax comp computed maximum principal stress in speci- 
men j based on P^ and assumed material 
parameters 

o 0 scale parameter 

a u threshold stress (location parameter) 

00 characteristic strength 

Subscripts 

assumed assumed 

comp computed 


1 



previous 

j 

T 

v 

s 


assumed previous value of assumed 

specimen designation 

tensile 

volume 

surface 


Analysis Based on Three Parameter 
Uniaxial Weibull Model 


The three parameter uniaxial Weibull model is used to 
describe the material inert strength probabilistic distribution. 
For both volume and surface flawed specimens, least squares 
best fit (LSBF) methods are developed to obtain the three 
material parameters from experimental tests on pure tension, 
pure bending, and four or three point loaded specimens. The 
necessary and sufficient condition for a solution is satisfied 
when the three computed parameters produce the lowest 
value of the sum of the residuals squared, that is, when 
n \2 

X IGfi - (jfj ] = minimum, where n is the number 

\ Umax. comp Umax } 

j=l 

of specimens tested. The (Pfj.Ofj^ ) data points are obtained 


from the experimental tests where Pq = (j - 0.3)/(n + 0.4). 
Pjj and^fj^ are, respectively, the probability of failure and 

maximum principal tensile stress in the specimen at fail- 
ure. Ofi is the computed maximum failure stress based 

Umax .comp 1 

on the value of Pfj and the computed inert strength material 
parameters. 


Pure Tension (Fig. 1(a)), Volume Flaws 


Pfj =1- exp 


G fj (x,y,z)-g uv 


\m v 


dV 


( 1 ) 


Gf](x,y,z) > g uv 

where Vjj is the volume in tension of the j ^ specimen with a 
stress distribution throughout the volume denoted by 
G^(x,y,z), G uv is the threshold stress, G ov is the scale factor, 
and m v is the Weibull modulus. For this case Gf](x,y,z) = 
a Qmax ^d tens ^ e S a ge v °l ume °f specimen j is Vjj = 
L^bjWj. Hence 

= m v ln(Gfj mjw -G uv )-m v lnG ov (2) 


In 


'n(»- p f ])'' 


V Tj 


The following system of n linear equations is solved in a 
LSBF sense: 


In 


In 


In(l-Pfi) 


-1 


*T1 


ln(l-P f2 ) 


-1 


V T2 


In 


InQ-Pfar 1 

V T „ 


f ln (°n raajt - CT uv) l 
ln (°f2 m « -^uv) 1 


'"K™ -°uv) 1 


m v 

m v lnG ov 


(3) 


In matrix notation { Y } = [A] {X}, where the j lh term in the 

ln(l - Pfj) -1 

column vector {Y} is yj = In 


(X} = 


m v 

[~m v lnG ov 
obtain the LSBF solution is 


1- 


and vector 


The equation that must be satisfied to 


{X} = [a t a] _1 [a t ]{Y} 


(4) 


where superscript T defines the transpose. 

The answer is obtained in the following manner: Assume a 
value for g uv » and solve for m v and g ov . With these values, 
compute the model failure stresses a fj maxcorap at all of the (n) 
Pfj data points, where 


Gfj ^ov 

U max, comp w 


ln{l — Pq) - ' 


+ G„ 


(5) 


Evaluate the sum of the squares of the residuals, where 

n 2 

Sum = V [Of; - a f| ) (6) 

\ Umax .comp Umax ) 

j-1 

Repeat the process for another value of G uv . Compute the 
new sum of the squares of the residuals (eq. (6)). Continue 
until the parameters (m^Go^o^) produce the minimum 
value of the sum of the residuals squared. 
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Pure Tension (Fig. 1(a)), Surface Flaws 


Pf) = l-exp 


i 


r (T lj (x,y,z)-CT u 
i \ G os 


dA 


where 8^ = 0 U v w j/( 2a Q max ) and V-jj = L 2 bjWj/2. The 
following expression is derived from equation (10): 


In 


(7) 


a Q (x,y,z)>a us 

For pure tension, Gfj(x,y,z) = °fj max , and the area in tension 
is Ajj = 2 L 2 (Wj + bj) where L 2 is the gage length, W is the 
width, and b is the thickness. Hence, 

ln(l-P rj ) 


In 




'Ti 


+ lnCT fw =0 + m v)m(a fw -ct uv ) 


+ In 


(l + m v )a 


m v 

ov 


( 11 ) 


In 


\-l 


A X j 


m s ln ( a f] max °os) m s lna os (®) 


Equation (8) is the basis of a LSBF evaluation of the Weibull 
parameters. From the inert data, a set of n linear equations is 
obtained. In matrix notation { Y } = [A] {X} where the j th 
term of the column vector {Y} is In [In (1 - Pq^VAtj] and 

vector {X}is< s l. The matrix [A] is the same as 
[-m s In G os J 

that in equation (3), except that the subscript v is replaced by 
s. The solution is obtained by the same method as that for the 
volume flaw solution. Assume o us , and solve for m s and . 
With these values, compute the failure stresses, Or, 

r max, comp 

for all n specimens. The computed failure stress for the j th 
specimen is 


The j* 1 term of the column vector {Y} is In [In (1 - Pfj) -1 / 

[ ( 1 + m v) 

Vrj] + In o f ; and the vector { X } is < v m i-i | 

J Jm “ [ln[(l + m v )ff“ v ] 

The matrix [A] is the same as in equation (3). Assume o uv , 
and from equation (11), solve the system of n linear 
equations in a LSBF sense, where j varies from 1 to n. In 
matrix notation, { Y } = [A] {X}, and the solution to this set 
of equations is obtained from equation (4). The final solution 
is the set of parameters associated with the minimum sum of 
the squares of the residuals defined by equation (6). They are 
obtained by the following procedure. Assume o uv , and solve 
for m v and o ov . With these values compute the predicted 
failure stresses <^Q max at all (n) Pfj data points. A simple 
method is to assume Or, . and solve for Or, 

■Umax .assumed ‘J max .comp 

where 


0f Jn 


= o„ 


1 

cT 

l 

N 

m s 

+ ® us 

(9) 


(l + m v )ln(l-P rj )'‘ 


max, comp 

A T j 



1 

> P 

1 


l + m v 




Evaluate equation (6), the sum of the residuals squared. 
Continue the process for another value of o us . Compute the 
sum of the residuals squared. Continue until the parameters 
(m^c^Ous) produce the minimum value of the sum by 
equation (6). 

Pure Bending (Fig. 1(b)), Volume Flaws 

Substituting the expressions dV = L 2 bj dy and Gg (x,y,z) 
= 2 tffj raax y/Wj into equation (1) results in 


(!-Pfj) ' = exp 


W: 


f \ ^ 




b; 


’ i2 l 


20fj y 

max 


v W i 

6n v J 


= exp 




/ \l+m v 

V Tj ( a fjmax _a uv) 


} o\ J 


1 + m. 




dy 


( 10 ) 


( 12 ) 

The next assumed value is Or, = 0.5 o^ + 

max. assumed I ‘J max, comp 

Or, . Repeat this process until Or, is 

*J max, previous assumed J U max. comp 

within some specified tolerance of Or, . Then 

compute the sum of the residuals squared by means of 
equation (6). Repeat all of the previous steps until the 
minimum value of the sum of the residuals squared is 
obtained. The parameters (m v ,Ob v ,c uv ) associated with this 
minimum are the solution. 

Pure Bending (Fig. 1(b)), Surface Flaws 

From equation ( 1 ) 


(l- p fj) 1 =exp 


1 ^ mS 

— f A (CT fj (x,y,z)-CT us ) ms dAj 

o rtv j “At: 


(13) 
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where Ajj is the tensile surface area of specimen j. 
Therefore, considering both the side and bottom surfaces of 
the specimen yields 



— m s ln(ofj m ^ a us ) l n [^os s (l + m s)] 


Assume a value for Of: . and iterate until o fl is 

Umax .assumed Umax, comp 

within some specified tolerance of Of; . Obtain the 

r U max. assumed 

sum of the residuals squared by equation (6). Repeat the 
process. The parameters that produce the minimum value of 
the sum of the residuals squared are the solution. 

Four Point Bend Specimen Fig. 1(c)), Volume Flaws 


Substituting the inner span and outer span stress distribu- 
tions a fj (x.y.z) = 20 f][nM yAV j and ^^(x.y.z) = 
4 °fj maJl Xy/(LiWj) into equation (1 ) results in 



where 5ij = OuvWj/^CJfj^ ). From equation (17) we obtain 


Solve in a LSBF sense the set of linear equations obtained by 
means of equation (15) and denoted in matrix form by 
{ Y } = [A] {X}. The j th term of the column vector { Y } is In 
[In (1 - PfjH/O^WjM - In [(1 - cj^/Of jmtt ) + (1 + m s )b/Wj] 


m s 

. f m /. \1 y The LSBF par am- 

ln[a“ s (l + m s )]j y 

eters are obtained in the following way: Assume a us , and 
keep this value Fixed. To evaluate the vector {Y}, assume a 
value for m s (m s assume( i) based on the two parameter 
solution. Evaluate the column vector { Y } and the matrix [A]. 
Solve for the vector {X} by equation (4). Compare the 
computed value of the Weibull modulus rn s?comp with 
n*s, assumed * Repeat this process until both values, m s comp and 
m St assumed» are within some specified tolerance. With these 
parameters (m^o^a^) and Pfj, compute the n values of 
Or, where 

v max, comp 

1 




(16) 





(18) 


Solve the set of linear equations obtained by equation (18), 
denoted in matrix form by {Y} = [A] {X}. For constant 
values Li and L 2 , the j* term of column vector { Y } is 
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The j th row of matrix [A] is [inGg^ l.ojand {X} = 
m v 1 

i T m v . \1 } . Assume a value for the threshold 

ln[<J ov ( 1 + m vjJJ 

stress a uv and an appropriate value for the Weibull modulus 
m v,assumed (based on the two parameter solution). Evaluate 
column vector {Y} and matrix [A]. Solve for solution vector 
{X} by equation (4). With C5u V fixed, solve for m vcomp . 
Iterate until m v cornp is within a given tolerance of m Viassumec ( . 
To compute the sum of the squares of the residuals by 
equation (6), we obtain °fj maxcorop from the following 

equation: 


Three Point Bend (Fig. 1(d)), Volume Flaws 

Substituting Ofj mM (x,y,z) = 4a fw /('L,W J ) into equation 
(1) results in 


« Lhi 
2 2 


in(l - Pfj) * = 2 bj J J 


r 4q fw x y 

L.Wj 


S|j 


-a„ 


dx dy (21) 


where 8^ = o uv Wj/(2o^ max ) and 8^ = L{WjC uy /(4Cf }m ^y). 
With V-jj = L]bjWj/2, integration of equation (21) results in 



°fj 


max, comp 


Vt-: 


l/m v 


(l + m v )ln(l — Pf] ) 


-1 


Wj /2 ( 

I! 


\l + m v 


8 >j 


2y 

W j aq 


max, assumed 


dy + 


L 2 


\l+m v 




max, assumed J 


(19) 


Assume a value of Or. , and iterate until is 

U max. assumed ‘Jmax, comp 

within some specified tolerance of o f . . Evaluate the 

Umax .assumed 

sum of the residuals squared by equation (6). Repeat the 
process. The parameters that produce the minimum value of 
the sum of the residuals squared are the solution. 

If all failures occur within the inner span and the tensile 
stress distribution outside the inner span is neglected 
(Lj = 0.0), equation (17) becomes the pure bend solution, that 
is, 


— 1 — exp 


L 2 bjWj 


2(l + m v ) 


Of. 

Umax 


\m v f \l +m v 

l__£uv_ 


V °ov J 






Umax / 


Therefore, 


In 


ln(l-P fj )- 

L 2 bjWj 


+ lnCI Qmaji -( 1 + m *) ln K aB " q uv) 

-ln[(l + m v )cC] (20) 


Equation (20) is the same as equation (11). 


In 


Infl-Pof 1 


W J 


V Tj 


-In 




5 U 


1 2y 

ylwj 0i] 


\l + m v 


dy 


max J 


= m v lna fjraM - ln[(l + m ] 

( 22 ) 


Equation (22) is the limit case of equation (18) with L 2 = 0.0. 
Solve the set of linear equations denoted in matrix notation 
by { Y} = [A] {X} by the LSBF method (eq. (4)). The jth 
value of column vector { Y} is 


yj = In 


4-n.r 


w j 


■In 


'Tj 


I 2 y G 

y l Wj 


\l+m v 


dy 


max / 


and the vector is {X} = 




-.To obtain an 


initial value of m v , let a uv = 0 and solve for the uniaxial 
Weibull two parameter distribution satisfying equation (4). 
Starting with this computed value of m v as m v assume j and a 
fixed value of a uv , evaluate the integral (eq. (22)) in column 
vector {Y}. The integral is evaluated numerically between 
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the lower limit 8jj and upper limit W j/2. Obtain from solution 
vector {X}, rn vcomp . A solution is obtained when m v comp is 
within some specified tolerance of m v assumed- When this 
does not occur, the next choice for m v assurne( j is 0.5 (m v comp 
+ m v»previous assumed)’ Iterate until m v>assume d is within some 
specified tolerance of m v comp . To compute the sum of the 
residuals squared by equation (6), we evaluate a fj rawcorop in 

the following way: For the n data values (Pq^fj^ ) where 
j = 1. ", assume Of WassunBd = Of inmi - The lower integration 
limit is 5jj = CT uv Wj/(20 5iMxtssuined ). With this limit, solve 
for a fi 

'J max.comp 




^OV^Ti 
max .comp A J 



(l + m 

V )ln(l-P fj )" 

5. 



2 

( 

\ 

f 1 

2y 

<*UV 

Jy 

«ii 

l W J 

G (j max .assumed j 


Assume a new value forOfj 


max, assumed 


= 0.5 ( Ofj 


l + m v 

dy 

(23) 

max .previous assumed 

+ CTfi ). Repeat the process, integrating over the new 

limit S xj until the previous assumed value is within some 
specified tolerance of the computed new value ( a fj maxcorap ~ 

Or, ). Compute the sum of the residuals squared by 

U max .assumed 

equation (6). Repeat the process assuming a new value for 
o uv , and continue until the minimum sum of the residuals 
squared by equation (6) is obtained. When this occurs, the 
values of m v , o uv , and o ov are the three material parameters. 

Four Point Bend Specimen (Fig. 1(c)), Surface Flaws 

Substitute into equation (7) the inner span side surface 
uniaxial tensile stress distribution Oq (x,y,z) = 20^^ y/Wj , 

the bottom surface tensile stress distribution o^(x,y,z) = 
Or, and the outer span uniaxial tensile surface stress distri- 
butions Of](x,y,z) = 4 Of )m ^ Xy/(LjWj) and a^(x,y,z) = 
2 Of, x/Li. Normalizing the area with respect to LjWj 

Umax J 

results in the following equation: 


Pg = 1 - expi 


x dy + 


W; 



p 

L,Wj 

til 

f 

2y 

A 

G us 

V a os 7 


(l + m s) 

J y 

5 

v w i 

a 0_ ; 


V+m. 


L 2 W j +L 1 b j 


L,Wj 


j _ ^us 


Nl+m s 


0 fj„ 


(l + m s )L 2 bj 


V Umax ) 


L| Wj 


1 _ G us 


V 




max / 


(24) 


where 6|j =WjCT us /(2 ). Thus 


In 


mj'-Pf,)’ 1 

- In 

"j 

2 ri| 

^y 

<>US 

LiWj 

J y 

8 -j ' 




\ l+m s 


+ 


L 2 Wj +L,bj 


LiWj 


dy 


o llc 1 S (l + m s )L 2 bj 


V 0f Jr„« } 


L,Wj 


< a 

J _ ^ us 

V °Cmax / 


= m s ln °5 M1 " ln °os 5 0 + m s )] (25) 


Solve the set of linear equations obtained from equation (25) 
(denoted by (Y) = [A] {X}, in matrix form) by the LSBF 
method (eq. (4)). The j A value of column vector { Y} is 


yj 


; = In 


l„(l-p,) _l 

-In 

TT j 

Ul 

t 

2y 

\ 

<*US 

LiWj 

J y 

V 

, w i 

a fjma* , 


1+m, 


L 2 Wj + Ljbj 


\1+m. 


L,Wj 


| _ g US 


V ^nuu ) 


dy 


(l + m s )L 2 bj 

LiWj 


( A m - 

] °us 

V af Jmax J 
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and vector {X} = 



+ m s 


)] 


j.Row j of matrix [A] is 


t a Qmax * ’0], Assume a value for the threshold stress a m and 
an appropriate value for the Weibull modulus, m s assumec j 
(based on the two parameter solution). Evaluate column 
vector {Y} and matrix [A]. The vector {X} is evaluated by 
equation (4). With a us fixed, solve for m s comp . When 
ro^comp i s within a given tolerance of m s assumec j, a solution 
results. To compute the sum of the squares of the residuals by 
equation (6), the values of a fl are obtained from the 

J max .comp 

following equation: 


with 5)j = OusWj/aaq^ ). 83 y = Lj Wj/(4o fw ), and 
6 3j - a us LyZOfj^.Thus 


In 


In 


LjWj 


-1 


- In 


W: 


z 

J 


2y q us 

V W j ° f Jraax J 


l+m t 


bj 

dy + — — 

W: 


\l+m. 


a fi™ J 


m s Ina fj m » - ln [( 1 + m s )a 


a fj 


max,comp 



(28) 


l 


m s 


(l + m s )ln(l-P g ) ' 


W j 


2 if 


v 1 + IDt 


1 


, 5 u 


V 


2y_ 

Wj a f] 


dy + 


L 2 W J +L 1 b j 


max, assumed J 


LjWj 


\l + m s / 


"A 


max, assumed 


LjWj 




max, assumed 


(26) 


To determine Of; , the process is the same as that 
outlined for the four point bend, volume flawed specimen. 
Likewise, the evaluation of the material parameters is the 
same as that outlined for the four point bend, volume flawed 
specimen. 

Three Point Bend Specimen (Fig. 1(d)), Surface Flaws 


In matrix form { Y} = [A] {X} and the jh term of column 
vector { Y } is 






'"('-Po) -1 

— In 

\ 1 

2y 

LjWj 

1 

1 



\l+m. 




Umax J 


dy 


From equation ( 1 ) 


ln(l - Pfj) - ' = 


( , \ m ‘ 


V °OS J 


Wj Lj 


J } 

5 u § 2,(y) L,Wl 




dx dy 




} ( 2xc fj ^ ms 

+ 2bj J — ^-<*us 


L; 

5 3j v J 


dx 


(27) 


+ 


W: 


\l + m c 


O n 


CT c 


Umax / 


and {X} = 



, \1 ?■ The solution to this set of 

+ m s)jj 


linear equations denoted symbolically by { Y} = [A] {X} is 
solved by equation (4). Starting with = 0.0, solve for m s 
and the scale factor Next, assume a value for The 
integrand is a function of m s . Starting with an assumed value 
of m s , iterate until m s assuraec t is within some specified limit of 
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comp. Then evaluate the scale factor. To find the value of 
<Jfi associated with the probability of failure Pf; , and 

1 J max, comp 1 J 

computed material parameters (m s ,a^a us ), satisfy the 
following condition: 


Omax,c< 




2 

/ \ 

f 1 

2y ^us 

J y 

S.j 

Wj CTfj 

y J J max, comp y 


l + m. 


dy 


W: 




1+m,. 


(j max, comp J 
1 


^ 1 + m X 

V L i W J J 




(29) 


Assume values for CTf, and iterate until the left side is 

equal to the right side constant. Compute the sum of the 
residuals squared via equation (6). Repeat the process 
assuming a new value for a us , and continue until the mini- 
mum value of the sum of the residuals squared is obtained. 
When this occurs, the parameters (m s ,0^,0^) are the 
solution. 


Three Parameter Specimen Uniaxial 
Weibull Model 

This report deals with material properties that are 
independent of the component geometry. However, a simple 
model is often used to obtain the inert strength probabilistic 
distribution of a given component (refs. 2 and 3). The charac- 
teristic strength parameter Oq in this model is component 
dependent and is not a material property. For completeness, 
this model is briefly mentioned. This model equation for 
volume flaws is formulated as 


Experimental Applications 

The examples in this section make use of some of the 
developed equations. Inert failure data are analyzed from the 
following Modulus of Rupture (MOR) bar data: 

(1) Four point bend room temperature failure data of 
sintered silicon nitride, table I (ref. 1). All failures were due 
to volume flaws and occurred within the inner span. These 
data were also used for four point bend surface flaw analysis 
to illustrate the application of the developed equations. 

(2) Three point bend transverse annealed silicon carbide 
data at 1300 °C , table II (Private communication from Sung 
Choi and Jonathan Salem, NASA Lewis Research Center). 
All failures were caused by surface flaws. 

(3) Three point bend longitudinal annealed silicon carbide 
data at 1300 °C, table III (Private communication from Sung 
Choi and Jonathan Salem, NASA Lewis Research Center). 
All failures were caused by surface flaws. 

To develop confidence in the method developed in this 
report, comparisons were made with the pure bend results 
from reference 1 . The equations for the four point bend and 
three point bend specimens were then developed and 
programmed. 

Sintered Silicon Nitride (Pure Bend Analysis, Volume 
Flaws, Table I) 

Monolithic silicon nitride data (SNW-1000, GTE Wesco 
Division, table I) obtained from reference 1 are used to com- 
pare the results of various LSBF techniques. All of the data 
in table I contain failures that occurred within the inner span. 
In reference 1, pure bend loading (fig. 1(b)) was therefore 
assumed applicable to these data. The three material parame- 
ters were computed using Cooper’s method (ref. 4), a modi- 
fied LSBF approach (ref. 5), and the method developed 
herein. Table IV summarizes the results of the three tech- 
niques used in the analysis of these data and the results 
obtained herein of the two parameter model (a uv = 0). Fig- 
ure 2 is a plot of the data points and the cumulative Weibull 
two parameter distribution curve. Figure 3 is a plot of the 
data points and the cumulative three parameter Weibull 
distribution curve. 


Pfj = 1 " exp 



CT fj O uv ^ 

‘Jmax u 

m y 


l °0v ) 



(30) 


For surface flaws, subscript v is replaced by subscript s. 
Since the characteristic strength is not a material property, 
this model has its limitations. It is commonly used and is 
mentioned for completeness. 


Silicon Carbide (Three Point Bend Surface Flaws, 

Tables II and III) 

Table V summarizes the results obtained for the two and 
three parameter uniaxial Weibull models. The cumulative 
distribution curves for three point bend data (Private 
communication from Sung Choi and Jonathan Salem, NASA 
Lewis Research Center) for transverse and longitudinal 
annealed silicon carbide and the data points are plotted in 
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figures 4 to 7. The two parameter distribution curves are 
plotted in figures 4 and 6. The three parameter distribution 
curves are plotted in figures 5 and 7. 

Sintered Silicon Nitride (Four Point Bend Analysis, 
Volume Flaws, Table I) 

Table VI summarizes the results for the two and three 
parameter uniaxial Weibull models. Figure 8 is a plot of the 
two parameter cumulative distribution curve and data points. 
Figure 9 is a plot of the three parameter cumulative distri- 
bution curve and data points. 

Sintered Silicon Nitride (Four Point Bend Analysis, 
Surface Flaws, Table I) 

The four point bend volume flaw data are also used for a 
four point bend surface flaw analysis to illustrate the applic- 
ation of the developed equations. It is realized that these data 
are not representative of the physical problem. Table VII 
summarizes the results for the two and three parameter 
uniaxial Weibull models applied to these data. Figure 10 is a 
plot of the two parameter cumulative distribution curve and 
data points. Figure 1 1 is a plot of the three parameter 
cumulative distribution curve and data points. 


Discussion and Conclusions 

Solutions are obtained from inert failure data based on the 
minimizing of the sum of the residuals squared as a necessary 
and sufficient condition. There are programs to evaluate the 
three parameters fitted to the specimen uniaxial Weibull 
model (eq. (30), refs. 2 and 3). The characteristic strength, 
ae vi a parameter in this model, is not a material property but 

component dependent. The results obtained using this model 
are only applicable to that specific component made from the 
same test material. 

In this report, material property parameter estimation 
methods are developed based on the uniaxial Weibull model 
(eq. (1)). The parameters so obtained are applicable to any 
component made from the same test data material. For the 
sintered silicon nitride four point bend inert volume flaw 
failure data (table I), Cooper's method (ref. 4), a modified 
LSBF method (refs. 1 and 3), and the approach developed 
herein were used to minimize the sum of the squares of the 
residuals based on the pure bend solution (all failures 


occurred within the inner span). Comparing the results 
reveals that the largest variation of the sum of the squares of 
the residuals (table IV) was less than 3 percent. Further 
comparison of the results of the three methods indicated the 
approach used herein was slightly less conservative in the 
low probability of failure regions (Pf < 0.05) and slightly 
more conservative in the upper region (Pf > 0.25). Figures 2 
and 3 are plots of the data points and the computed 
cumulative distribution curves for the two and three 
parameter uniaxial Weibull models based on the pure bend 
solution. 

The results for the silicon carbide three point bend surface 
flaw data (tables II and III) from longitudinal and transverse 
annealed specimens are summarized in table V. Comparing 
the two and three parameter models reveals that there are 
large differences in the Weibull modulii and scale factors. 
The cumulative distribution curves are plotted in figures 4 
to 7. Superimposing the two and three parameter curves 
reveals small but significant differences because most 
designs are based on the very low probability of failure 
region. 

The four point bend solutions to the two and three 
parameter uniaxial Weibull models applied to the data in 
table I are summarized in table VI. The cumulative distribu- 
tion curves are plotted in figures 8 and 9. Figures 2 and 3 are 
the cumulative distribution curves for the pure bend solution. 
A comparison of figure 2 with figure 8 and figure 3 with 
figure 9 indicates the four point bend results for both cases 
are slightly more conservative in the lower probability of 
failure region and slightly less conservative in the higher 
probability of failure region. 

Four point bend volume flaw data in table I are used for a 
four point bend, surface flaw analysis to illustrate the 
application of the developed equations. It is realized that 
these data are not representative of the physical problem. The 
results are plotted in figures 10 and 11 and summarized in 
table VH. 

Justification for applying the three parameter model rather 
than the two parameter model will depend on which model 
better predicts the behavior of a component with its observed 
performance. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, March 1996 
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TABLE I — FOUR POINT 
BEND SILICON NITRIDE 
VOLUME FLAW INERT 
FAILURE DATA (fig. 1(a)) 
[Li = 20.8 mm, L 2 = 19.6 mm. 


TABLE II. — THREE POINT BEND SILICON 
CARBIDE SURFACE FLAW TRANSVERSE 
ANNEALED INERT FAILURE DATA 
(fig. 1(c)) 

[Span = L 1 = 19.936 mm.] 




Specimen 

number 

Thickness, 

b } , 

mm 

Depth, 

Wj. 

mm 

Failure 

load, 

kg 

Specimen 

number 

Failure strength, 
MPa 



613.9 


1 

2.991 


m 


623.4 


2 

2.999 

1.875 

H 


639.3 


3 

2.999 

1.873 

16.20 


642.1 


4 

2.999 

1.877 

14.85 

5 

653.8 


5 

3.000 

1.874 

15.08 

6 

662.4 


6 

2.995 

1.871 

13.13 

7 

669.5 


7 

2.996 

1.875 

14.18 

8 

672.8 


8 

2.998 

1.874 

15.00 

9 

681.3 


9 

2.997 

1.879 

14.78 

10 

682.0 


10 

2.999 

1.875 

11.63 

11 

699.0 


11 

2.997 

1.874 

12.15 

12 

714.5 


12 

2.998 

1.877 

14.33 

13 

717.4 


13 

2.998 

1.877 

14.33 

14 

725.5 


14 

2.999 

1.879 

13.23 

15 

741.6 


15 

2.995 

1.876 

12.83 

16 

744.9 


16 

2.994 

1.874 

15.75 

17 

751.0 


17 

2.998 


16.23 

18 

761.7 


18 

2.997 

1.876 

12.83 

19 

763.9 


19 

2.994 

1.877 

12.75 

20 

774.2 


20 

2.997 

1.876 

13.05 

21 

791.6 


21 

2.999 

1.879 

16.05 

22 

795.2 


22 

3.001 

1.875 

14.85 

23 

829.8 


23 

3.000 

1.877 

16.23 

24 

838.4 


24 

3.000 

1.875 

13220 

25 

856.4 


25 

2.993 

1.871 

17.63 

26 

868.3 


26 

2.993 

1.879 

12.30 

27 

882.9 


27 

2.995 

1.879 

16.05 



28 

2.996 

1.876 

15.08 



29 

2.994 

1.877 

13.05 



30 

2.996 

1.877 

10.05 



31 

2.996 

1.872 

12.75 



32 

2.995 

1.874 

17.70 



33 

2.996 

1.876 

15.30 



34 

2.994 

1.878 

11.55 


TABLE III.— THREE POINT BEND SILICON 
CARBIDE SURFACE FLAW LONGITUDI- 
NAL ANNEALED INERT 
FAILURE DATA (fig. 1(c)) 

[ Span, L 1 = 19.936 mm. 1 


Specimen 

number 

Thickness, 

bj, 

mm 

Depth, 

Wj, 

mm 

Failure 

load, 

kg 


2.999 

1.866 

wgm 


2.998 

1.866 

R9 


3.001 

1.867 

17.18 


2.998 

1.868 

15.53 

5 

2.998 

1.869 

14.10 

6 

3.000 

1.863 

15.68 

7 

3.000 

1.871 

15.98 

8 

3.003 

1.870 

16.95 

9 

3.001 

1.863 

17.25 

10 

3.002 

1.863 

12.00 

11 

2.991 

1.866 

15.75 

12 

2.994 

1.864 

15.23 

13 

2.992 

1.866 

14.78 

14 

2.993 

1.866 

14.25 

15 

2.995 

1.868 

14.93 

16 

2.996 

1.869 

12.15 

17 

2.9% 

1.870 

15.38 

18 

2.9% 

1.871 

14.78 

19 

2.997 

1.871 

14.70 

20 

2.997 

1.872 

12.90 

21 

3.002 

1.873 

9.98 

22 

3.001 

1.874 

10.95 

23 

3.000 

1.867 

15.53 

24 

3.002 

1.871 

14.40 

25 

3.003 

1.864 

13.50 

26 

3.000 

1.864 

12.53 

27 

3.000 

1.864 

12.90 

28 

3.000 

1.865 

13.95 

29 

3.002 

1.864 

14.70 

30 

3.002 

1.863 

12.45 

31 

2.992 

1.859 

12.75 

32 

2.993 

1.868 

13.95 

33 

2.992 

1.868 

15.53 

34 

2.994 

1.869 

17.55 

35 

2.975 

1.870 

12.08 
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TABLE IV.— WEIBULL PARAMETERS OBTAINED FROM SILICON NITRIDE (SNW-1000) FOUR POINT BEND VOLUME 

FLAW INERT FAILURE DATA (TABLE I) 

[Data in table I are analyzed as a pure bend solution over the inner span.] 


The probability fo failure for a given value of cr fj max i s defined as Pfj = 1.0 - exp 


v Tj "l(°fjmax ~ q uv) 


l + m v 


V 1 + m v 


°f)n 


Pure bend solution 
(volume flaws, fig. 1(b)) 

Weibull 

modulus, 

m v 

Scale factor, 

Oov* 

MPa-m 3/m v 

Threshold stress, 

duv, 

MPa 

Sum of residuals squared, 
27 - 

Yk - o t ) 

\ computed data / • 

j = l 

Starlinger, et al. and Cooper - LSBF (ref. 1) 

1.625 

0.00258276 

560.84 

2684.4 

Starlinger, et al. - Modified LSBF (ref. 1) 

1.677 

0.00370464 

558.08 

2664.0 

LSBF a 

1.608 

0.00218469 

565.195 

2743.7 

Two parameter LSBF methodb 

11.306 

150.1733 

0.0 

13440 


a LSBF applied to eq. (11). 
b a uv set equal to zero in eq. ( 1 1). 


TABLE V.— WEIBULL PARAMETERS OBTAINED FROM THREE POINT BEND SILICON CARBIDE 
SURFACE FLAW INERT FAILURE DATA (TABLES II AND III) 


The probability of failure for a given value of Ofj^ defined as 


Pfj = 1.0 - exp j- 


°fJn 


V CT os 


(l + m s ) 


J - 

8 ,j y 


\l+m. 


2y 


V W J CTf Jmax; 


dy + ' 


W: 


\ 1+ nic 


1 - - 


°5 


max / 


where 6jj = 


^usWj 


2a n 


Umax 


LSBF best method 
(surface flaws, 
fig. 1(d)) 

Weibull 

modulus, 

ms 

Scale factor, 
MPa - m 2 /m & 

Threshold 

stress, 

Ous, 

MPa 

Sum of residuals 
squared, 

i (of -o f ) 

^ computed data J j 

Transverse- annealed two- 
parameter model 

9.294 

114.52 

0.0 

2665 

Trans verse-annealed three- 
parameter model 

4.024 

11.71 

190.0 

1942 

Longitudinal -annealed two- 
parameter model 

9.161 

114.14 

0.0 

1417 

Longitudinal- annealed three- 
parameter model 

5.893 

39.78 

120.0 

1259 



TABLE VI. — WEIBULL PARAMETERS OBTAINED FROM SILICON NITRIDE (SNW-1000) FOUR 
POINT BEND VOLUME FLAW INERT FAILURE DATA (TABLE I) 


The probability of failure for a given value of G q is defined as 

TW: 


P f =l-exp- 


r L l b jWj 


CT ov J 2(l+m v ) 


J ~\ 

6 lj y 


2 y q uv 

w i a fi 
J r J max j 


l+m v 


dy+ 


z2t\ 


N l+m v 


1— 


^uv 


°fjmax 


where 8 jj = 


q uv w j 


2 °fjr 


LSBF method 
(volume flaws, 

fig. 1(c)) 

Weibull 

modulus, 

m v 

Scale factor, 
Oov, 

MPa - m 3/m v 

Threshold stress, 

0|1V» 

MPa 

Sum of residuals 
squared, 

21 i \2 

\ 0i J max.comp °Qmax ): 

j = l' 'J 

Weibull two 
parameter model 

10.841 

141.713 

0.0 

11712.5 

Weibull three 
parameter model 

1.443 

0.0006804 

564.0 

2038.5 


TABLE VII.— WEIBULL PARAMETERS OBTAINED FROM SILICON NITRIDE (SNW-1000) FOUR POINT BEND 

VOLUME FLAW INERT FAILURE DATA 

[For illustrative purposes the data in table I are analyzed as surface flaw inert failure data.l 


The probability of failure for a given value of Of. is defined as 

max 


Pf. = 1 - exp* 


( Of. ^ 

Omax 

^ ^os ) 


LiWj 


(l + m s ) 


f - 

V 


( 


\ 1 + m<j 


2y_ 

w j °fj 


max J 


dy + 


L 2 Wj + Ljbj 


( 


LjWj 


\l + m s 




max J 


u ^ a us W j 

where 5^ = — 


2o (jr 


(l + m s ) 


L 2 bj 


\">s 


L,Wj 


°fj 


max J 


LSBF method 
(surface flaws, 
fig. 1(c)) 

Weibull 

modulus, 

m* 

Scale factor, 
aos, 

MPa - m 2/m* 

Threshold stress, 

Ous, 

MPa 

Sum of residuals squared, 

21 i \2 

^ l^Qmiu,comp ^Qraax )■ 

j = l J 

Weibull two 
parameter model 

10.84 

325.23 

0.0 

12067.0 

Weibull three 
parameter model 

2.00 

1.69 

575.0 

1928.3 





Fj/2 F/2 

(d) 


Figure 1 . — Specimen loading and geometry, (a) Pure tension, (b) Pure bend silicon nitride (SNW-1000) specimen; 
|_2 = 19.6 mm, b = 4.0 mm, W = 3.1 mm. (c) Four point bend silicon nitride (SNW-1 000) specimen; Li = 20.8 mm, 
L 2 = 19.6 mm, b = 4.0 mm, W = 3.1 mm (table I). (d) Three point bend silicon carbide specimen; L| = 1 9.936 mm 
(tables II and ffl contain Fj, bj, and Wj values). 
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Failure stress, <r f , MPa 


Figure 2. — Distribution curve for two parameter Weibull model from volume flawed inert siiicon 
nitride (SNW-1000) data via specimens subjected to pure bend analysis over inner span (Fig. 1(b), 
table I). Weibull modulus m v = 1 1 .306; scale factor <r ov = 150.173 MPa-m 3/mv ; threshold stress 
tr uv = 0.0 MPa; sum of residuals squared = 13440. 



Failure stress, <r f , MPa 

Figure 3. — Distribution curve for three parameter Weibull model from volume flawed inert silicon 
nitride (SNW-1000) data via specimens subjected to pure bend analysis over inner span (Fig. 1(b), 
table I). Weibull modulus m v = 1 .608; scale factor <x ov = 0.0021847 MPa-m 3/mv ; threshold stress 
cr uv = 565.2 MPa; sum of residuals squared = 2744. 
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Figure 4. — Distribution curve for two parameter Weibull model from surface flawed inert transverse 
annealed silicon carbide data via specimens subjected to three point bend analysis (Fig. 1(d), 
table II). Weibull modulus m s = 9.294; scale factor cr os = 114.5 MPa-m 2 ^™ 3 ; threshold stress 
tr us = 0.0 MPa; sum of residuals squared = 2665. 



Failure stress, o fl MPa 


Figure 5.— Distribution curve for three parameter Weibull model from surface flawed inert transverse 
annealed silicon carbide data via specimens subjected to three point bend analysis (Fig. 1(d), 
table II). Weibull modulus m s = 4.024; scale factor a os = 11.708 MPa-m 2 *" 8 ; threshold stress 
a us = 190.0 MPa; sum of residuals squared = 1942. 
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Figure 6. — Distribution curve for two parameter Weibull model from surface flawed inert longi- 
tudinal annealed silicon carbide data via specimens subjected to three point bend analysis 
(Fig. 1(d), table TII). Weibull modulus m s = 9.161; scale factor a oa = 114.14 MPa-m^ 018 ; 
threshold stress cr us = 0.0 MPa; sum of residuals squared = 1417. 



Figure 7. — Distribution curve for three parameter Weibull model from surface flawed inert longi- 
tudinal annealed silicon carbide data via specimens subjected to three point bend analysis 
(Fig. 1(d), table III). Weibull modulus m s = 5.893; scale factor a QS = 39.78 MPa-m 2 ^ 3 ; 
threshold stress o us = 120.0 MPa; sum of residuals squared = 1259. 
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Failure stress, a f , MPa 


Figure 8. — Distribution curve for two parameter Weibull model from volume flawed inert silicon 
nitride (SNW-1000) data via specimens subjected to four point bend analysis (Fig. 1 (c), table I). 
Weibull modulus m v = 10.84; scale factor a uv = 141 .7 MPa-m 3/mv ; threshold stress 
ct uv = 0.0 MPa; sum of residuals squared = 11713. 



Failure stress, a f , MPa 

Figure 9. — Distribution curve for three parameter Weibull model from volume flawed inert silicon 
nitride (SNW-1000) data via specimens subjected to four point bend analysis (Fig. 1(c), table I). 
Weibull modulus m v = 1.443; scale factor <r ov = 0.000680 MPa-m^ mv ; threshold stress 
c r uv = 564.0 MPa; sum of residuals squared = 2039. 
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Failure stress, tx f , MPa 


Figure 10. — Distribution curve for two parameter Weibull model from volume flawed inert silicon nitride (SNW- 
1000) data via specimens subjected to four point bend analysis. As an illustrative example these data were 
analyzed as if they came from surface flawed specimens (Fig. 1(c), table I). Weibull modulus m 8 = 10.84; 
scale factor tr os = 325.2 MPa-m 2 ^ 8 ; threshold stress a us = 0.0 MPa; sum of residuals squared = 12067. 



Failure stress, o f , MPa 

Figure 11 . — Distribution curve for three parameter Weibull model from volume flawed inert silicon nitride (SNW 
1000) data via specimens subjected to four point bend analysis. As an illustrative example these data were 
analyzed as if they came from surface flawed specimens (Fig. 1(c), table I). Weibull modulus m s = 1 .997; 
scale factor tr os = 1 .6895 MPa-m 27 ™ 8 ; threshold stress o us = 575.0 MPa; sum of residuals squared = 1928. 
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